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ABSTRACT 



UHECRs are roughly isotropic and attain very large energies, E ;> 3 x 10 20 eV. Conventional 
models fail to explain both facts. I show here that acceleration of UHECRs in GRBs satisfies 
both observational constraints. Using Meszaros and Rees' (1994) model of GRBs as due to 
hyperrelativistic shocks, I show that the highest energies that can be attained thusly are 
E ~ 10 20 9 ~ 5 / 3 n 1 5 ^ 6 eV, explaining the energy of the Bird et al. (1995) event even without 
beaming. The traditional photopion catastrophe affecting UHECR acceleration in AGNs is 
circumvented. An order of magnitude estimate shows that the total energy flux of UHECRs 
at the Earth is also correctly reproduced. A test of the model based upon the UHECRs' 
distribution upon the plane of the sky is briefly discussed. 



Subject headings: acceleration of particles - gamma-rays: bursts 
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1. Introduction 

Ultra High Energy Cosmic Rays (UHECRs, E > 3 x 10 18 eV, the 'ankle') show a flatter spectrum 
than Cosmic Rays (CRs) just below the ankle. They have energies extending up to an observed maximum 
of E max = 3 x 10 20 eV (Bird et ai, 1995), are thought to be predominantly protons (Bird et al, 1994) 
and arrive from directions roughly isotropic on the plane of the sky. Since expected angular deflections for 
highly energetic protons are less than 10° (Sigl, Schramm, Bhattacharjee 1994, SSB from now on), this 
rough isotropy is thought to reflect a (rough) intrinsic isotropy of their sites of production. 

The conventional acceleration mechanism invoked so far is first order Fermi (1949) acceleration 
at shocks (Axford, Lee and Skadron 1977, Blandford and Ostrikcr 1978, Bell 1978, Krymsky 1977). 
Acceleration sites must be located within the finite range (<; 100 Mpc, SSB) over which UHECRs can 
travel before their energy is significantly degraded by energy losses due to photopion and photoclcctron 
production off CMBR photons. However, all acceleration sites proposed so far (sec SSB for a beautiful 
review) fail on two counts. First, they cannot attain the highest energies observed so far, falling short by at 
least one order of magnitude. Second, they cannot reproduce the rough isotropy of the directions of arrival, 
since they explain all UHECRs as coming from the handful of peculiar objects (AGNs, radiogalaxies, the 
Virgo Cluster and so on) that can be found within the finite range mentioned above. Even more frustrating 
is that no prospective candidate for an acceleration site can be located within a suitable error box around 
the direction of arrival of the E = 3 x 10 20 eV event (SSB, Elbert and Sommers 1995). 

We should thus be scouting around for a class of objects that is roughly isotropically distributed, and 
where sufficient amounts of concentrated energy are available to accelerate UHECRs. One such category, 
so far unexplored, is gamma-ray bursts (GRBs, see Paczyhski, 1993 for a review). I shall consider only 
cosmological models of GRBs, with a total energy release of Eqrb — 10 51 erg. I shall need in the following 
discussion no detailed property of the mechanism proposed to explain the energy injection mechanism, but 
I shall need the details of the hydrodymanical expansion of the fireball leading to the GRB. In particular, 
the really attractive feature that I shall try to exploit in the following is the suggestion by Meszaros and 
Rees (1994, MR from now on) that hyperrelativistic shocks (whether due to the impact of different parts of 
the same flow, endowed with different Lorenz bulk factors, or to the impact of the flow on the surrounding 
interstellar medium) are responsible for GRBs. 

In Section 2, I shall give two generic arguments (i.e., independent of the actual acceleration mechanism) 
in favor of GRBs as sites for the acceleration of UHECRs. Those features of the hydrodynamics of GRBs 
which are relevant to the problem at hand, are briefly reviewed in Section 3, which is entirely based upon 
the results of Meszaros and Rees (1994, and references therein). In section 4, I describe the acceleration 
process in this model for the evolution of the fireball. In particular, in Section 4.1 I discuss qualitatively 
two acceleration mechanisms, and I compute in subsections 4.2 and 4.3 the highest energies that UHECRs 
can attain. A mixed bag of limitations and caveats are discussed in Sections 5. I discuss the results in 
Section 6, and summarize them in Section 7. 



2. Gamma ray bursts as accelerators of UHECRs 

There are two arguments that make GRBs appealing, the first one being a numerological coincidence. 
The total energy of UHECRs striking the Earth can be estimated as 5 x 10~ 13 erg s _1 cto~ 2 sr^ 1 . This 
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has been obtained by taking the UHECRs' spectrum as N(E) cx E~ 2 - 7 , with the normalization coming 
from Hillas (1984). This local flux must be compared with that released in the form of UHECRs by GRBs. 
Since UHECRs have a finite range (< 100 Mpc, SSB), I consider only the nearby GRBs, which occur at 
the rate of hp = 30 yr~ 1 Gpc~ 3 (Paczyhski 1993), each releasing about Eijhecr — 10 51 erg in the form of 
UHECRs. The total energy striking the Earth (if the sources are uniformly distributed, and are all standard 
candles) is thus hpD m EijHECRl^ = 4 x 10~ 13 erg s _1 cmT 2 sr , independent of beaming. Here I took 
D m = 100 Mpc as the maximum range for E = 3 x 10 18 eV, the threshold energy above which I computed 
the observed total energy flux of UHECRs striking the Earth. I have assumed a kind of cquipartition, such 
that the energy released by the GRB in the form of photons equals that in the form of UHECRs, probably a 
not too bad assumption, given the large relative velocities freely available in GRBs. It also agrees with the 
well-known high efficiency of particle acceleration at strong shocks (Volk, Drury and McKenzie 1984). The 
above coincidence of theoretical and observed fluxes is quite striking. It is, to the best of my knowledge, the 
first time that the amplitude of the UHECRs' flux at Earth has been 'explained'. 

The second argument that makes fireballs attractive is that they can act as hyperrelativistic ping-pong 
bats with respect to CRs. GRBs are known to be highly super-Eddington: in fact, they show substructure 
in their pulse profile on a scale ~ 1 ms, implying source sizes < 3 x 10 7 cm, which are, at most, the 
Schwarzschild radii of 100 M Q objects, for which the Eddington luminosity is ~ 10 40 erg s _1 . Since GRB 
luminosities are 10 51 erg s _1 if at cosmological distances for typical burst durations of 1 s, clearly they are 
super Eddington. It seems likely thus that radiation pressure drives a relativistic expansion, with large 
Lorenz factors 7&. This conclusion is strengthened by the following argument. GRB spectra extend at 
least to 100 MeV. Even more astounding is the recent discovery (Hurley et al., 1994) that at least one 
of these spectra extends to 18 GeV . If this energy derived from thermal energy, then electron/positron 
pair creation would easily make the burst optically thick, and the emerging spectrum would be thermal 
(Paczyhski, 1986). Since instead GRB spectra are known to be much broader than black-body, 7& must 
refer to bulk motion kinetic energy (Paczyhski, 1993). The hydrodynamic evolution of a fireball shows that 
a shell containing all of the fireball's energy and baryon content is accelerated to relativistic Lorenz factors 
7i,, and remains thin through most of its evolution. A magnetic field B, initially at cquipartition, decreases 
as the shell expands; when the fireball ejecta impact upon the surrounding interstellar medium, the shocked 
shell magnetic field is revived to equipartition values. Thus, the expanding shell contains a significant, 
random field B which might deflect incoming cosmic rays backward. We would then have a very efficient, 
first order Fermi acceleration, for consider a highly relativistic particle of Lorenz factor -fcR moving radially 
inward toward the origin, which, after having entered the shell, may be deflected backwards by a magnetic 
irregularity which is comoving with the expanding flow. In the reference frame of the flow, it can easily be 
shown that relativistic composition of velocities implies that the cosmic ray has Lorenz factor ~ ^blCR 
both before and after being turned backwards, but after the deflection the cosmic ray, now moving radially 
outwards, has, with respect to the lab frame, a Lorenz factor ~ ^il^icR'- i-e-, the single backward deflection 
has increased its energy by a factor 47 2 ~ 10 5 , for the values jb — 10 2 — 10 3 favored by MR. Thus, if a few 
such cycles can be achieved, it seems possible that CRs can be accelerated up to energies of 3 x 10 20 eV, 
the highest energy event detected by the Fly's Eye. 

These generic remarks, which are meant to be independent of the specific acceleration mechanism, 
make GRBs palatable. Below I shall discuss the physical processes involved in the deflection/acceleration 
mechanism. 
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3. The hydrodynamics of GRBs 



I review here briefly some properties of GRBs which are relevant to the acceleration of UHECRs. This 
subsection is entirely based upon the results of Meszaros, Laguna and Rees (1993, MLR from now on). Rees 
and Meszaros (1992) have proposed that thermalization of the spectrum of the GRB can be avoided by 
introducing a slight contamination of baryons within the original fireball, parametrized by the parameter 77 

where M is the total baryon mass. MLR showed that, for a certain range of 77, only a fraction of the initial 
energy is still in the form of radiation at the moment in which the expanding fireball becomes optically thin, 
the rest having been converted into kinetic energy of the baryons. This kinetic energy becomes available for 
radiation in the optically thin regime if the baryons' directed kinetic energy can be suitably randomized, 
either by collisions with the external interstellar medium, or with slower or faster portions of the relativistic 
flow itself (Meszaros and Rees 1993). Since it seems rather difficult and contrived to produce a fireball with 
very little baryon contamination, and since this is the favoured model of MR, below I shall concentrate 
exclusively on high-load fireballs, i.e., those with 

1< V < T m = 3.3 x W 5 El{ 3 r- 2/3 , (2) 

where -E51 is the total energy release Eq in the fireball in units of 10 51 erg, and r§ is the radius ro in which 
such energy is released in units of 10 6 cm. For these values of 77, MLR showed that the baryons can be 
accelerated up to 7& = 77, ending up with just about all the initial energy E in kinetic form, while only a 
fraction r)/T m is released in photons in the first, miniburst. After an initial period (which is of no interest 
to us) in which the fireball accelerates from rest up to 7b = 77, a period of free expansion follows. In the lab 
frame, the energy in this phase is concentrated in a thin slab, initially of roughly constant thickness dr ~ ro, 
equal to the radius of the region in which the initial energy deposition took place. The free expansion 
begins when the shell has reached in the lab frame a radius r ~ r s such that 

r s = 77<r Vo (3) 

where 9 is the beaming semi-opening angle, 0—1 corresponding to isotropic emission, and ends when, in 
the lab frame, the shell starts expanding, linearly with radius r. This occurs for r ;> r&, with 

n = V 2r o ■ (4) 

I shall need in the next section the expression for several quantities in the frame comoving with the shell, 
expressed in terms of the lab frame shell distance from the origin, r. MLR give, for the comoving radiative 
energy and temperature 

JL\ - (L\ - J T 1/3 ^ 2/3 (H)A-) 2/3 r s <r<r b 

Eo) \To) \ V 1/3 0- 2/3 (ro/r) r > r b W 

Here T = 4 x 10 11 K El( 4 r 6 3 ^ 4 is the initial temperature. They also give the shell's thickness in the 
comoving frame dr as 

5r = { ™ ^< r < r " (6 ) 
{ r/rj r > r h 
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MLR also argue that a large class of phenomena leads to the growth of a magnetic field B, in near 
equipartition with the initial radiative energy. They deduce that the comoving value of B decreases 
according to 

B__( 7T 2/3 0~ 4/3 (n)A) 4/3 r s <r< r b 

and the initial, equipartition field Bq is given by 

B a = 1 17 El( 2 r 6 3/2 e /2 G (8) 

where £ is the ratio of magnetic to radiative energy, and measures the departure of B from equipartition 
(which occurs for £ = 1). 

The release of the energy stored as baryons' kinetic energy, in the form of the 7-ray flash occurs either 
by collisions of the parts of the fireball's endowed with different values of the Lorenz factor (Meszaros and 
Rees f993), or by impact of the fireball ejecta on the interstellar medium. The two different mechanisms 
have been incorporated into a single, coherent picture (MR) to explain the exceptional event observed by 
Hurley et ai, 1994. I shall consider only the second, latter shock, and shall speak loosely of shocked shells 
to refer to this scenario only. The physics of these collisions (Meszaros and Rees 1993, MLR) is similar to 
that of the development of SN shocks. A collisionless shock moves into the unshocked material with Lorenz 
factor V2jb, while a reverse compression wave moves into the shell to be decelerated, eventually steepening 
into a mild shock with a modest Lorenz factor, ~ 2. The material immediately behind the forward shock 
is heated to a thermal Lorenz factor as 7&. In the shocked shell, several phenomena (MLR) can revive the 
magnetic field to equipartition. 

The impact of the fireball ejecta occurs at a deceleration radius 

r d = 10 18 O-^E^n^rj- 2 / 3 cm (9) 

where the particle density per cm 3 in the ISM is n\. The shocked shell has thickness 

r sh = 10 18 0" 2 / 3 i4 /3 nr 1/ V 5/3 cm (10) 

and I scale the intershell magnetic field B with the value reached by equipartition with the baryons' energy 
density, which is the same, whether the material has been shocked by the forward or reverse shock. I find 

B = {).bn- 1,2 V 6- 1 t 1 l 2 G . (11) 

where £ again parametrizes departures from equipartition (£ = 1). MLR also consider the possibility that 
B does not, after all reach equipartition. In that case, it is limited from below by the fossil magnetic field 
(Eq. 0) formed at the outset of the expansion. For simplicity, they show that this corresponds to setting 
£ w 4 x lCT 3 in Eq. |ll|. 



4. The acceleration of UHECRs 



4.1. The acceleration mechanism 
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There are actually two distinct acceleration mechanisms that I want to consider. The first, possibly 
less interesting one is as follows. The freely expanding shell can deflect backward some CRs which, after the 
scattering, as discussed in Section 2, shall be boosted up in energy by a factor 4^ ~ 10 5 . In this mechanism, 
the shell simply scoops up whatever CRs are floating around, acting rather passively. Clearly, the flux 
estimate of Section 2 does not apply to this mechanism. However, it is nonetheless interesting because, 
after the explosion that shall eventually lead to the GRB, since the shell expansion is hyperrelativistic, 
business outside goes on as usual, and it is known by direct observations that Cygnus X-3 (Samorski and 
Stamm 1983), the Crab pulsar (Dzikowski et ai, 1981, 1983), Hercules X-l and Vela X-l (see Protheroc 
1994 for references and a critical, but optimistic review) produce cosmic rays with E « 10 16 eV. The 
presence of high-energy cosmic rays around normal neutron stars also agrees with theoretical prejudice, 
since most conventional estimates (Hillas 1984, SSB) of the highest energy attained by cosmic rays around 
neutron stars agree on the value E ~ 10 17 eV. Thus, if any of these CRs were to move backward and land 
on the expanding shell, it would be boosted up to very high energies. This mechanism can be extended to 
a self-consistent one, by considering the possibility that suprathermal particles (possibly preexisting the 
fireball explosion) are accelerated by repeatedly scattering off two subportions of the fireball, each having 
given rise to a relativistic shell. It is easy to see that the highest energy thusly attainable is the same as in 
the scoop-up model. The details of the scattering mechanism limit the highest energy attainable: this is 
discussed in Section 4.2. 

The second, more appealing acceleration mechanism occurs because the MR model for GRBs explicitly 
predicts the existence of hyperrelativistic shocks, because of the impact of the expanding flow on the 
surrounding interstellar medium. Since such events are seen to have a complex structure, it seems likely 
that in any case we are witnessing the collision of several subportions of the flow with each other, so that the 
GRB consists of several sections of converging flows, alternatively separated by a reverse shock, a contact 
discontinuity, and a forward shock. The contact discontinuity is likely to be unstable, leading to mixing of 
material on either side of it, so that particles from the extreme Boltzmann tail of Lorenz factor J> 7^ are 
injected over the whole shocked shell. The suprathermal particles shuffle between the pre-shock shell and 
the post-shock shell, each time being scattered by the shell's magnetic field (Eq. Q). At every such loop 
their energy is boosted up by a factor 47^ w 10 5 (see the argument of Section 2 and, most importantly, 
the simulations of Quenby and Lieu 1989). Even if we take the energy of suprathermal protons to be as 
low as 76 (i.e., equal to that of thermal particles), we see that after two cycles the protons have energy 
167^ > 10 , roughly enough to account for the highest energy event observed so far. Thus, the proposed 
mechanism is essentially a two-cycle Fermi-Bell acceleration. This mechanism is of course limited by the 
ability of the shell to deflect UHECRs backward. Limitations arise because the shocked shell has both 
finite thickness and finite lifetime, both to be considered in the following Section 4.3, together with other 
sundry problems. It is however clear at this point already that the highest energy that can be attained 
is reached when the shocked shell is largest, and the whole shocked shell thickness is equal to one mean 
free path for the UHECR's deflection. At the same time, another subportion of the flow has collided with 
this shell, and it too is completely shocked. Then the CR shuffles back and forth between two shells of 



thickness given by Eq. 10 each with a revived magnetic field (Eq. |llj). The computation below shall use 
this argument. Also, I should add that this is a coherent mechanism, one i.e. where the injection and 
acceleration of suprathermal particles are all due to the same physical environment, so that the estimates 
of the efficiency of the generation of cosmic rays (Volk, Drury and McKenzie 1984) and the equipartition 
argument of Section 2 apply. 



I now compute the maximum energy to which CRs can be accelerated. Such upper limit exists because 
the shells have a finite thickness and magnetic field B so that an UHECR's mean free path to diffusion, 
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being a multiple of the particle's gyroradius, cannot exceed the shell's thickness; outside the shell the 
magnetic field is in fact negligible with respect to that inside the shell. The computation shall be carried 
out in the shell's reference frame, and then the limiting energy shall be transformed to the lab frame. A few 
important comments on the physics of the scattering of CRs are deferred to the Section 5, in order not to 
obstruct the flow of the argument. 



4.2. The interaction of CRs with the expanding fireball 

I begin with determining the optical depth to photopion destruction of an incoming CR, as a function 
of the shell's expansion radius. Using ov = 1CP 28 cm 2 as the relevant cross-section (Caldwell et ai, 1978), 
I find (with the computation done in the comoving frame) 



E 

t v = n y a n dr = — 

and, using Eq. [g, I find 



T - = n ^ = kTAxr** *** (12) 



TV = 1.3 x 10 14 (1°) 2$V /4 . (13) 
Optical thinness to photopion destruction is then achieved beyond r„, 

rv = 1.1 x 10 7 r E 3 5 { 8 r~ 5/8 . (14) 
This occurs for r < or for r > rj, depending upon whether rj > rji, or 77 < rji respectively, where 

m = 3.3 x lO 3 ^ 716 ^ 5/16 . (15) 

Once the shell has become optically thin to photopion destruction, every incoming relativistic CR finds 
itself in a magnetic field, with respect to which it is super-Alfvenic. In fact, in the comoving frame where 
baryons are locally at rest, the Alfven speed Va, given by V\ = B 2 /4irp, is determined by the baryons' 
rest-mass density pi,. This is because, as shown by MLR, for r > r s , the internal energy in the shell 
comoving frame is subrelativistic (this is really the meaning of the parameter r s ), and the same applies for 
an equipartition B-field. Thus any CR with jcr J> 2 is clearly super-Alfvenic. Explicitly, we have 

V_l = _B^_ JE = J , ? 2/3 - 2 /3 r < n 

c 2 AirpbC 2 Mc 2 1 E Q ^ \ ^/a^/sro r>n \ > 

where I used Eq. [|. The CR being super-Alfvenic, the excitation of the usual helical modes that lead to 
scattering is possible before r ~ r n : in fact, Va/c < 1 for r > r s . In particular at the time of thinning to 
photopion destruction, r/r is given by Eq. pTl, so that I obtain 



UJ r m > v > m 

For low-load fireballs, 77 > T m , I derive the maximum Alfven speed at photopion thinning, 



(17) 



1 /3 

^ = O.Og^ 1716 ^ 5748 ^ 72 f — ] = 0AE^ 7/Ui r^ /72 e /2 (18) 



c \ rn 



- 9 - 



independent of 77, which shows the shell to be subalfvenic for any load. 

The deflection mean free path of relativistic cosmic rays is generally taken, by analogy with the 
interplanetary medium, to be a multiple g of the particle's gyroradius r^, where g ~ 40 for relativistic shocks 
(Quenby and Lieu 1989). Thus, in the comoving frame, the deflection of a CR with energy E = 10 15 £?i5 e V 
requires that the shell's thickness (Eq. |6|) is at least g times its gyroradius, i.e., gr L < Sr. Using 

rL = lpC W^G)' (19) 

and Eq. [?] and |] we get an upper limit to the comoving energy of the CR which can be deflected backwards. 
I find 



,. , , - 1.0 x lOV 773 ^ 7 2 r- 1/2 _4/3 I ( ^ )4/3 (2°) 

10 15^1/2 e y ' 51 6 \ r± r>n 



This gives the highest CR energy that can be deflected backwards by the shell after the thinning radius 
r n . Before that time, a smaller range of CRs could still be deflected, provided the total optical depth to 
photopion destruction seen along their walk through the shell, ~ grL, is less than unity; in other words, 
provided T n grL/Sr < 1. This, and Eq. |l3|, ^|, |l^ give 



lQl5£l/2 e y 



6.4 x 10-V /3 £5i l7 % 3/ V 4 / 3 { fe) ' r<n (21) 

r > r b 



From Eq. |2fj and |2l] it can be seen that the highest energy in the comoving shell is attained just at the 
moment of optical thinning r T , and, in the lab frame, this maximum energy is given by 

4/3 



1/4-1/12, -4/3,1/2 I {m) 7]>m 
51 '6 C \ f x 2/3 



E max ~ 2riE max {r„) = 4 x 10 iy eV E^r^'^e-^^ { ) "( 2/3 (22) 



This peak on the maximum energy for r = r^ is due to the competition between photopion destruction, 
which dominates at small radii, and the decrease in the comoving magnetic field, which makes the CR's 
gyroradius increase beyond the shell's thickness. 

Another loss mechanism that is potentially important is due to synchrotron radiation, which damps the 
UHECR energy on a timescale t s — 1 yr (10 20 eV/E)(l G/B) 2 . This is to be compared with the time that 
the UHECR spends within the shell, where the field is high, t& = Igr^jc. Using the formulas above, I find 

k_ui^(i)' (23) 

for the highest energy CR that can be reflected at shell radius r. At r = rv, we have t s /td ^> 1, and 
even more so later on, which shows synchrotron losses to be negligible, when no beaming is present. 
Alternatively, the synchrotron cooling time provides a lower limit to the amount of beaming consistent with 
the acceleration of UHECRs, since we must have t s /td > 1. Supposing 9 < 1, there is a characteristic radius 
rg such that synchrotron losses do not damp the CR, given by 

r» = 2^!. (24) 
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If rg < rv, the highest energy attainable is given by Eq. £2], otherwise, going back to Eq. |2(j, I find 

, f CJL W 3 n > 716 

E {e) -2nE (r a ) - 9 x 10 19 eV P l / 2 E l/2 r~ l/2 I Kr >e > , 1 ^ e 1 ' 2 (25) 

^max — ^'l^maxV 8) — » * lu e f S -^51 '6 S / _J)_ ^2/3/1-1/3 ^ t?o 



where I have conveniently defined t\q = (0.3^ 3//8 ) 1 / 2 77; . The above equation is the result we were searching 
for. It shows that UHECRs are easily produced in the free expansion phase of the events leading to a GRB, 
for a moderate beaming. In fact, using MR's favoured value r\ ~ 10 3 , I find that a beaming of ~ 0.01 ~ 1° 
is necessary to produce the highest energies observed to date, ~ 3 x 10 20 eV. 

It is interestng to notice that, since -fb < r m ~ 4 x 10 5 in GRBs, whether high or low load, (Shemi 
and Piran 1990, MLR), there is an absolute, universal maximum to the CR energy in this acceleration 
mechanism, given by 

E abs = 1.9 x 10 22 < /9 r 6 - 5/ V 4 / 3 ^ 2 eV . (26) 



4.3. The acceleration of UHECRs in the shocked shell 



The acceleration of UHECRs in the shocked shell after the impact of the fireball ejecta on the 
interstellar medium is subject to the same physics as in the previous paragraph. The Alfven speed is, as 
remarked at the end of Section 3.1, below the speed of light and the photopion catastrophe is no longer a 
problem. In fact, the energy radiated is now the whole kinetic energy of the baryons, so that the optical 
depth to photopion destruction is 

TV = = 1 X (™*) 2#V/» f^L_) 4/3 ( 27 ) 

e4wr 2 d \ e J 51 1 1000/ 1 ; 

where e is the average photon energy in the lab frame. If all of the energy were released by synchrotron 
losses, e would be in the X-ray (MLR), giving tv ~ 0.1. However, most of the losses occur by Inverse 
Compton (MLR) in the 7-ray band, giving t,<1. 

Proceeding as in the previous subsection, and equating g times the Larmor radius to the shell thickness 
I find the maximum energy in the comoving frame 

^(com.sh) 

10 r 5 a ; ve /2 = 5000 *- b/ v a/8 i$ v /a (28) 

and, in the lab frame, 

E^l = 2 V ElZT sh) = 10 19 O^'W^W" eV . (29) 

Synchrotron cooling is not a limiting factor, leading to a very weak upper limit on nj. I find here 

^ = 900 nfV 3 / 3 i4 /3 C 3/2 • (30) 
td 

In order to account for the Bird et ai, 1995, event, we must have 

0-5/3^1/3^1/3^5/6^1/2 > 3Q ^ 
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which gives 

and, since it is necessary that t s > td 



^ < 0.1 n- 5/ V 3 / 15 i4 3/15 r 1/5 (32) 
td 



n, < 0.08r? 26 / 25 i4 6/25 r 1/5 cm" 3 . (33) 
For r\ ~ 10 3 , the above limit becomes rt\ <J 100 cm -3 , which is a very weak limit. 
For the MR's favoured value r\ ~ 10 3 I find 

E^l = l^El{\-^9-^e' 2 eV (34) 

which has a sufficiently steep dependence on 8 and n\ to accommodate the energy of the 3 x 10 20 eV event 
without any trouble, even in the case of no beaming. In particular, it seems obvious that energies up to 
~ 10 22 eV are quite compatible with reasonable values of m. This is so especially because at least one 
(Narayan, Paczyhski and Piran 1992) of the mechanisms proposed to explain the initial release of energy 
predicts such events to occur outside galaxies, and all (Usov 1992 and 1994, Thompson and Duncan 1993) 
use pulsars, which are not confined to the Galactic disk: n\ is accordingly reduced. Then values as low as 
£ i=s 10~ 3 are compatible with the energy of the Bird et al. (1995) event. Such low value of £ also allows the 
possibility that the shell's magnetic field necessary to accelerate UHECRs is provided completely by the 
fossil magnetic field (Eq. |?j): it shall be remembered from the end of Section 3 that this field corresponds 
to £ « 4 x 10~ 3 . 

Eq. |34| is the major result of the paper. 

Emax is subject to an absolute upper limit since rj < T m , just like the limit for the freely-expanding 
fireball, derived in the previous section. I find 

= 7 x lO 20 ^! 7 V 7 V 5/ V /2 0~ 5/3 eV . (35) 



Lastly, it should be mentioned that the shocked shell has a finite lifetime, but that this provides no 
limitation on the maximum energy. In fact, the acceleration process has, as a bottleneck, the time that the 
UHECR spends on its last trip before being scattered for the last time, because its mean free path increases 
with energy. Thus, the acceleration occurs (in the shell frame) on a timcscalc sa 5?"l/ c = Sr/c which equals 
the light shell-crossing time. This is of course the shortest timescale on which the GRB can be generated, 
and thus the shock lasts at least as long as this (MLR). 



5. Caveats 



The main emphasis of this paper is on computing the highest energies that can be attained by UHECRs 
in GRBs, and in fending off the most obvious loss mechanisms, synchrotron and photopion. This explains 
the cavalier treatment reserved to relativistic shock acceleration of UHECRs. Below I try to make this 
treatment plausible. 

Acceleration of CRs at relativistic shocks has been studied by several authors, both in the test particle 
regime (Peacock 1981, Kirk and Schneider 1987, Quenby and Lieu 1989) and in the nonlinear regime (Bell 
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1987, Jones and Ellison 1991), and even in oblique shocks (Kirk and Heavens 1989, Ballard and Heavens 
1991). In the above, I basically used a test particle approximation, which could fail if the modification of the 
shock structure due to the inclusion of CRs were such as to decrease the relative speed of the incoming and 
outgoing streams. However, simulations by Bell (1987) and Jones and Ellison (1991) clearly display 'thin' 
shocks even in the nonlinear regime, and Bell (1987) shows that the detailed shock structure is irrelevant 
for the highest rigidities. 

In the previous section, I have assumed that the cosmic ray velocity is subject only to pitch-angle 
scattering, and that it is reversed in each scattering, which may seem unrealistic. Quenby and Lieu (1989) 
argue, on the basis of the analogy with the interplanetary medium, that scattering is essentially isotropic, 
and trajectory integrations for the fully relativistic case (Moussas, Quenby, and Valdes-Galicia 1987, 
Valdes-Galicia, Moussas and Quenby 1992) have shown that scattering occurs through large pitch-angle 
changes, with <5(cos6*) ~ 0.5 — 1. From this they deduce, through their numerical simulations, that the 
energy is increased by a factor 7^ per cycle, when the proper average over all cosmic rays' velocity directions 
is taken. In relativistic shocks the diffusion approximation breaks down because the cosmic rays' velocities 
are not isotropically distributed, but are instead strongly peaked toward the radial direction (see Fig. 1 of 
Quenby and Lieu 1989). Since I too assumed a radially peaked velocity distribution, the energy increase 
7^ applies to my case as well, and the difference with my previous computation, 47^, is all due to the 
substitution of backward-forward scattering with isotropic (although large pitch-angle) scattering. This 
argument thus validates my use of forward-backward scattering, the only ensuing error being a modest 
factor 4 in the energy increase per scattering. 

The special relativistic Alfven speed is given by 

Yl _ B 2 /4n 

~J = 4e 6 /3 + £ 2 /47r (36) 

where e;, is the (relativistic) baryons' energy density, and the factor 4/3 becomes 1 in the non-relativistic 
limit. Scaling to equipartition values I find Va/c w £ 1 / 2 /2. I argued above (see discussion after Eq. 
that even £ w 10~ 3 is acceptable, yielding Va/c » 0.01. However, for equipartition values the Alfven speed 
is w c/2, more than in all simulations mentioned above. Unquestionably, in this limit I am stretching 
the applicability of the usual turbulence arguments to the boundary of the non-relativistic regime. The 
most critical limitation arises in the assumed, phenomcnological link between the mean free path to CR 
scattering, A, and its gyroradius, A ~ gr^ = 40r^ (Quenby and Lieu 1989). Still, one should keep in mind 
that, in the model of the previous section, matter shocked by the reverse shock, in the shell's reference 
frame, is barely relativistic, roughly as assumed by Quenby and Lieu (1989), and that relativistic effects for 
Va/c~ 0.5 do not appear so extreme to force one to abandon the previous estimate (Eq. p4[ ). 

Next, one may wonder whether sufficient strong turbulence is present to ensure scattering of UHECRs 
in the shocked shell. I have two arguments about this. First, it seems quite likely that, given the large 
velocities and energies available, very strong magnetic turbulence can be generated behind the shocks. 
Second, this magnetic turbulence may have been observed already, albeit indirectly. It is well-known, in 
fact, that GRBs are extremely rich in substructure down to a scale of w 1 ras. In MR's model, since the 
GRB arises from synchro-Compton radiation, this substructure may be interpreted as inhomogeneities of 
the magnetic field: where the magnetic field is strongest, the generation of synchrotron radiation is more 
effective, and its conversion into 7-ray radiation leads to the local peaks in the observed time structure 
of the GRB. Whether this substructure corresponds to small-scale shocks, instabilities, or nonlinear wave 
effects, it must correspond to inhomogeneities in the magnetic field on the corresponding wavelengths. As a 
matter of fact, this argument is similar to that usually made (Quenby and Lieu 1989) to justify the presence 
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of large-scale inhomogeneities in the magnetic field strength in the jets of radiogalaxies. There, however, 
synchrotron radiation emission is directly observable, while in GRBs it must be inferred from the 7-ray 
observations. 

The magnetic field associated with the shell is of course impossible to predict theoretically. The only 
guide we can invoke is the interstellar medium, where the magnetic field is well-known to have achieved 
equipartition. On the other hand, the initial mass motions must be so violent, that any combination of 
compression, shearing, turbulent dynamo, magnetic buoyancy (Usov 1992 and 1994, Narayan, Paczyhski, 
Piran 1992, Thompson and Duncan 1993) are likely to lead to near equipartition values. The same 
comments apply of course to the shell after it has started to decelerate, so that the achievement of 
approximate equipartition seems to me by far the most reasonable assumption. It should be noted, however, 
that MLR also consider the more extreme case of 'magnetic dominance', in which the magnetic energy is in 
equipartition not just with the energy to be released eventually in the GRB, but with the binding energy of 
the object leading to the fireball. In this case the frozen-in magnetic field would a factor of 30 higher than 
in Eq. |^, in which case the estimeates of the maximum energy (Eq. ^5|) would be an underestimate by the 
same factor of 30. 



6. Discussion 

I discuss briefly here some consequences of the hypothesis that UHECRs are connected with GRBs. 

Can we be sure that most UHECRs are protons, rather than nuclei with higher charges, as suggested 
by observational evidence (Bird et at, 1994)? If GRBs are extragalactic, the gyroradii of iron nuclei in the 
intergalactic field (~ 10 _9 G, SSB), are of order IMpc, for a 10 20 eV particle. Thus, they are confined to 
large bubbles surrounding their sites of production, and rather far from us. Then, diffusion in the Galactic 
magnetic field, in which their gyroradius is 1 kpc for the same energy, tends to bar them access to the inner 
regions of the Galaxy where we are observing them, and carry them outwards. Thus, if GRBs accelerate 
UHECRs, the problem of their composition is automatically solved, even neglecting photodestruction of 
heavy nuclei at the acceleration sites. 

One should not expect a close temporal association between UHECRs and GRBs. Naively, one might 
think that, since the UHECRs' speed differs from c by one part in Jq R <; 10 19 , UHECRs should trail GRBs 
by less than 10 -3 s, even if they arrive from about 100 Mpc away. This would lead to the expectation 
that, every time we see an UHECR, 7-ray satellites should observe a GRB. However, UHECRs are bent 
along their path by the intergalactic magnetic field by approximately 10° (SSB), leading to a path longer 
by ~ 10 15 s, thus washing away any correlation with GRBs. 

The expected angular distribution of UHECRs is isotropic by construction (because such are the 
GRBs, Meegan et at, 1992), and an absolutely unavoidable consequence of the model. Hopefully, it ought 
to become testable with the construction of the Giant Airshower Detectors. Minor departures from exact 
isotropy are however expected. If GRBs are cosmologically distributed, the expected dipole is of order 
1 — 2 x 10~ 2 , but most importantly it is oriented in the direction of the Sun's peculiar motion with respect 
to the frame of the CMBR (l,b) = 264°. 7, 48°. 2 (Maoz 1994). However, in the cosmological case the 
expected dipole for UHECRs would depart from that of GRBs because energy losses prevent the arrival of 
UHECRs from distances ;> 100 Mpc (SSB). This is especially interesting because it is exactly in this region 
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(Scaramella, Vettolani and Zamorani 1994) that the Local Group's peculiar velocity forms. In particular, 
this seems interesting especially in connection with the well-known anisotropy of UHECRs (Hillas 1984), 
which has often been ascribed to the contribution of the Virgo Cluster, which is also known to contribute 
greatly to the formation of the Local Group's peculiar velocity. Furthermore, UHECRs have a direct 
cosmological application in that, by restricting attention to higher and higher energy bins, we can select to 
look at closer and closer distances, and thus ought to be able to see the dipole moment fade away. These 
anisotropies are currently being studied and the results shall be published elsewhere. 

There is very little that can be said about the UHECRs' spectrum, at this stage. In the test particle 
case in relativistic shocks (Kirk and Schneider 1987), the CRs' spectrum is harder than in non-relativistic 
shocks. Non-linear effects can modify this conclusion, even though it is not clear in which direction: Bell 
(1987) argues that, if injection is limited to w 1CU 4 of the background number density, the test particle 
spectral shape is valid, while the non-linear computations of Jones and Ellison (1991) find a softer spectrum 
than for the test particle case. How each of these computation would be modified by a nearly relativistic 
Alfven speed is not known. This spectrum should also be convolved with the steepening induced by 
photopion and photoelectron destruction induced by propagation over distances J> 10 Mpc. 



7. Summary 

This paper is essentially about a coincidence. Conventional sites for the the production of UHECRs 
fail to explain both the highest energies observed so far, and the rough isotropy of their directions of arrival. 
It is suggestive that there is one class of objects, those that produce GRBs, which can remedy both these 
defects. 

The directions of arrival of UHECRs in this model are isotropic because such are GRBs (Meegan et al., 
1992). Two acceleration mechanisms are qualitatively described in Section 4.1. The highest energies that 
CRs can attain by bouncing off fireballs are given by E ~ 9 x 10 19 9-^ 3 eV when the fireball is in the phase 
of free expansion (see Eq. |25| ). When, however, we consider as an acceleration mechanism the conventional 
Fermi acceleration at shocks, and the shocks are the highly relativistic ones invoked by MR to explain 
GRBs, then the highest energies that CRs can attain, are given by E ~ 10 20 _5 / 3 ni 6 £ 1,/2 eV (Eq. |J). 
In the first case (no shocks), some beaming is necessary to explain the highest energies observed to date 
(Bird et al., 1995), while in the second one a proper choice of n\ is sufficient, thus doing without beaming. 
If the UHECRs are generated at the same time as GRBs, and if some kind of equipartition between the 
various forms of energy losses is achieved in nearby (<; 100 Mpc) GRBs, then the total flux of UHECRs at 
the Earth is correctly predicted by an order of magnitude estimate (see Section 2). Lastly, the photopion 
catastrophe has been shown to be irrelevant after the fireball has expanded beyond a radius r„ (Eq. |l4|) . I 
have argued that a test based upon the distribution of UHECRs' directions of arrival is feasible, and that 
no conclusion about the spectral shape is possible at the moment. 

Thanks for helpful discussions are due to R. Scaramella and E. Pesce, and especially to B. Paczyhski. 
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